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Abstract 

The effects of a strong magnetic field in the dynamical response of a pure neutron system to a 
weak neutrino probe are studied within the framework of the Landau Fermi Liquid Theory in the 
non-relativistic Hartree-Fock approximation. We use Skyrme forces to parametrize the partially 
magnetized nuclear interacting plasma and describe its collective modes. We study the vector and 
vector-axial matter response functions and the neutrino mean free path in this two-component 
fermion system with net spin polarization and density dependent in-medium correlations. We find 
a decrease in the neutrino opacity of magnetized matter with respect to the non-magnetized case 
for fields log w B(G) > 17. 
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I. INTRODUCTION 



The theoretical study of properties of hadronic matter under extreme conditions allows 

n 

to understand nuclear systems ranging from experimental heavy ion collisions [1[ to the 
astrophysical scenarios with objects such as neutron stars, where densities up to several 
times that of nuclear saturation are believed to exist in their inner core The study 
of the low temperature behaviour of Fermi Liquids was pioneeringly addressed by Landau 
[3I and later developed by others in the Fermi Liquid Theory (FLT) 4|. The main idea 
underlaying is that the properties of this type of normal quantum systems can be studied 
using the dynamics of quasiparticles close to the Fermi surface. The inclusion of an additional 
magnetic field, B, allows further testing the properties of magnetized Fermi Liquids since 
nucleons have a tiny, but not vanishing, magnetic moment. 

These degenerate plasmas are of interest in the description of matter in the inner shells 
of rapidly rotating neutron stars with very intense magnetic fields. In the catalog of the 
~ 1,500 such objects there is now experimental indication that in a subgroup of pulsars, 
called magnetars, magnetic field strengths in the surface of the order B m 10 15 G 5), [f| 
exist. In the formation of these objects, after a Supernova core collapse, neutrinos are an 
important ingredient in the first stages of the cooling of the proto-neutron star The 
interactions of (anti) neutrinos with matter determine the ratio of deposited energy that, in 
turn, affect the dynamics of the cooling afterwards. 

As a first approach, and in most of the calculations, the modifications due to the induced 
spin polarization in nuclear systems under the presence of magnetic fields are discarded. For 
example, some works have calculated the matter response functions in absence of magnetic 
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is included 
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neutron system [£ , |9|, [lOj or in the case where some proton population 
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171 . Il8| . For systems where a spin excess is postulated 



■priori some partial results have been obtained in the context of Hartrcc-Fock calculations as 
in [l^]. However, in these works there is no link between a magnetic field and the induced 
matter magnetization. 

It has been argued that it is important to know whether a ferromagnetic transition takes 
place or not [2oj] in nuclear systems at high densities and the impact of the spin polarization 
as a source of additional excitation modes in the magnetized plasma that, in turn, could affect 
the neutrino opacity in the nuclear medium. To describe the nuclear interaction, a diversity 
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of models in the relativistic 



2l| or phenomenological non-relativistic [12J approximations 



have been considered. On the other hand, in the low density case, recent works show that 
the vector response seems to be relevant in the low density non-magnetized matter and 
direct simulations either using Monte Carlo or Molecular Dynamics techniques have been 



performed. At these low densities clustered pasta phases may develop [22], |23|, |24], 125 ]. 

In an earlier work we showed that there are important modifications in the low 
and high density regimes, as magnetization of a neutron plasma is concerned, when one 
considers phenomenological nuclear interaction models, such as Skyrme and Gogny forces. 
At supranuclear densities, 1—3 times saturation density, the residual magnetization amounts 
to values typically below 5% for magnetic fields as high as B m 10 18 G allowed by the scalar 
virial theorem and moderate temperatures. We have also calculated the influence of the 
generalized Landau Fermi liquid parameters [2^] on a set of observables in a magnetized 
neutron system [28 ]. 

The present work performs a study of dynamical response functions for a magnetized 
neutron plasma and the effect of a strong magnetic field on the neutrino opacities in the 
two-component neutron Fermi Liquid, using Skyrme effective nuclear interactions. The 
structure of the paper is as follows. In section [Til we describe the medium response under 
the linear approximation for the different parametrizations of the Skyrme force used in the 
description of the nucleon-nucleon (NN) interaction. Using the Landau Fermi Liquid Theory 
at zero temperature we consider a multipolar expansion of the particle-hole (ph) interaction 
and analyze the linear response up to dipolar terms. We gauge the effects of finite density, 
presence of a strong external magnetic field and the induced spin polarization in the neutron 
system. In section IHII we present the results obtained and finally summary and conclusions 
are given in section IIVI 



II. RESPONSE FUNCTIONS IN THE MAGNETIZED NEUTRON FERMI LIQ- 
UID 

We consider a homogeneous pure neutron system where each neutron has spin projection 
on the z-axis, a. It can be either a = +1 for spins aligned parallel to a uniform magnetic 
field that will be taken along the z-direction, B = Bk, or a = — 1 for antiparallel spins. 
This system may be viewed as a two-component neutron Fermi Liquid where the relative 
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populations of magnetized species vary as magnetic field and density change. Although we 
restrict ourselves to neutron matter in this work, for application to astrophysical scenarios, 
i. e., in neutron star matter [29j beta equilibrium must be imposed and, subsequently, 



multi-component plasmas with leptons, protons and heavier baryons should be considered. 

The total baryonic density in our two-component neutron system is given by p = p + + 
being p a = ^pr the cr-polarized component at T = 0. In this kind of systems there may 
be some spin excess, A, defined as A = p+ ~ p ~ . As explained in for given conditions of 



density, p, temperature (we fix here T = 0) and a strong magnetic field strength, B, the 
relevant thermodynamical potential is the Helmholtz free energy, Fm, defined as [301 ] 



M 



E-MB, 



(1) 



where E is the energy and M = J mdV is the total magnetization of a given volume. The net 
magnetization density is m = p n Ap and n n = — 1.9130427(5)/i7v is the neutron magnetic 



moment in units of the nuclear magneton 3lJ. Note that parallel (antiparallel) aligned 
magnetic moments (spins) are energetically favoured. 

In this work we have considered an effective approach to the nuclear interaction using two 
different parametrizations of the phenomenological zero-range Skyrme force. In the usual 
way it can be written as 32] 



V^ rme ( ri , r 2 ) = t (1 + XoP") S(t) + \t x (1 + x x P°) [k' 2 5(r) + 5(r)k 2 

+t 2 (1 + x 2 P CT ) k' • 5(r)k + \h (1 + x 3 P°) p a (R)6(r) 

o 



(2) 



where r = r x — r 2 and R = (ri + r 2 )/2, k = (Vi — V 2 )/2z the relative momentum acting 
on the right and k' its conjugate acting on the left. P a is the spin exchange operator. 
Note that we have omitted terms not relevant for homogeneous systems. From the large 
myriad of Skyrme parametrizations we have chosen the widely used SLy4 and SLy7, given 
by the Lyon group 
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34j as they provide good values for binding of nuclei and also for 



neutron matter equat 



around 1.5M P 
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ion of state (EOS) giving values of maximum neutron star masses 



361 ] . Saturation density is slightly different for both parametrizations, 
being for the SLy4 (SLy7) case 0.158 /m" 3 (0.160 fm~ 3 ). 



In the context of the FLT the properties of these systems at low temperature are related to 
the dynamics of quasiparticle excitations with momentum k and spin projection, a, around 
the Fermi surface in the presence of a magnetic field. These are studied by evaluating the 



quasiparticle interaction matrix element |4| and will be a crucial ingredient needed in order 
to compute the response of the pure andpossibly magnetized hadronic system. They can be 
obtained for a non-magnetized system 37| from the usual multipolar expansion in Legendre 
polinomials of the quasiparticle interaction and we can write, 

oo 

V ph = E \fi + gi^ 2 }Pi(cos9), (3) 

1=0 

where 9 is the angle related to the interacting quasiparticle three-momenta and and g% 
are the so-called Landau parameters. In the more general case where two possible spin 
orientations a = ±1 are taken into account, the polarized quasiparticle matrix elements can 



be written 



27j using coefficients j\ a ' a ^ depending on the degree of polarization, A, and on 



the Fermi momenta k F , k F involved. In the limit A — ► the Landau parameters fullfill the 
following relations, 

fi = , (4) 

Act, a) _ Act -a) 

9l = h f • (5) 

For the effective Skyrme interaction that we will be considering in this work the only non- 
vanishing terms are the monopolar (I = 0) and the dipolar (I = 1) ones. In the Skyrme model 
they are independent of temperature, so that they just retain the density and polarization 
dependence. They can be written for the monopolar case as (a' = ±cr), 

ft a) = lt 3 (l-x 3 )(a(a-l)p a - 2 p +P+ + 2ap a - 1 (p-p (7 ))- 
o 

+t 2 (l + x 2 )k 2 F>a , (6) 

ft' a) = t (l - z ) + k(l - x 3 )[a(a - l)p a ~ 2 P+P - + (a + l)p a ] - 

o 

" xi) + t 2 {l + x 2 )) (k 2 ^ + , (7) 

while the dipolar terms are given by, 

A (CT ' CT) = -t 2 (l + x 2 )4, CT , (8) 

ff'~ a) = ~(ti{l ~ xi) + Hi + x 2 ))k F ,M,-cr. (9) 

In the Skyrme interaction the quasiparticle effective mass in a magnetized system depends 
on the polarized dipolar coefficients, 



k 2 
I 



K/m = 1 + -NoAft' a) + (^ z£ )/ 1 (<J '" a) ] (10) 



where N 0a 



m*k 



F.tT 



2tt 2 



contains also m* and are the quasiparticle level densities at each po- 



larized Fermi surface. The response functions in a neutron system for a weakly interacting 
neutrino probe under the presence of a magnetic field B can be obtained from the Lan- 
dau FLT through the matter susceptibilities for a given multipolarity, I 38| . Accordingly, 



the Lindhard function, ^to°"') = x^ a ' a '\<^, g), satisfy the Bethe-Salpeter equation and can 
be written under the form of an algebraic system that, in the monopolar random phase 
approximation (RPA), reads as a 4 x 4 matrix, 



X 



to 



?{aa") <a"a') 



(11) 



where Xo — I (Cj3^o corresponds to the non-correlated response due to the a-polarized 
neutron fraction obtained from integration of the quasiparticle propagator Gq. In the limit 
where the energy transfer is close to zero, uj —>■ 0, and u/q\s approximately constant, then 
Re[xl\ -> -N 0a . 

The explicit solution in the RPA monopolar (I = 0) case can be written as, 

X i+ ' a ' } = T7 (x^^cr', +) - x^x^ [ft^S^, +) - ft^W, -)] ) , (12) 



X 



D 



D 



x^ 



(13) 



where the determinant of the set is 



(14) 



In the dipolar approximation (/ < 1) we can write the RPA algebraic coupled system in- 
cluding the susceptibility as a 8 x 8 matrix 23] 



X 



_ A?) 



x^O^O+x^ E ffV^ + iP E /^r^vo, (is) 



) x toV) +7 to E ^r^o, (i 6 ) 



where we have defined the following auxiliar quantities 



rto^) = y 



(27T) 



rCOS 



(0) G M , 



^=1^0.^8)^. 



(17) 
(18) 
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The expressions for the coefficients can be written in the Landau limit as 

l { i ] = v {a) x£\ (19) 

7 W = v^xP ~ (20) 

where v^ a > = ™°^ q - In-medium effects are incorporated by the polarized quasiparticle 
propagator, and corrections to that at the linear approximation through G^ a,r7 '\ ob- 

tained by summing all the ring diagrams corresponding to the magnetized ph excitations 
calculated at a given /-multipolarity. Then, the Lindhard function is calculated from the 
propagator by integrating over the incoming neutrino momentum as, 



^j^X (2!) 



(2tt)= 

Having obtained the generic reponse in the (a, a') channel, we can now constrain to the 
isovector contribution, S = 0, as, 

x (s=o) = + x (-) + + x (- +)j (22) 

and the vector axial contribution, 5=1, reads, 

x (s=i) = x ( + +) + x (--)_ x (+-)_ x (-+). (23) 

The structure functions can be calculated from the Lindhard function at zero temperature 



as 



38] 



S s=0 '\u, q) = —Im X 5=0,1 (w, q). (24) 

7T 

In this case of vanishing temperature only non-negative energy transfers may happen in the 
polarized neutron system. Then, the response function will provide information about 
the density-dependent in-medium correlations, partial polarization of the plasma, and vector 
and vector-axial contributions to the overall response of matter that may arise in the presence 
of a strong magnetic field when it is weakly probed by neutrinos. For example, in the neutral 
current reaction involving u^, v T neutrinos, collectively u x , and their antiparticles, elastically 
scattering up (down) polarized neutrons, 

v x + n a — > v x + n° . (25) 
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In the non-relativistic limit the differential cross section of neutrinos scattering off matter 



39] 



can be calculated from the linear response of the medium using 

' <i " Cf Gl -E'\l - f u {E')){C 2 {\ + cos9)S°(u, q) + C 2 A (3 - cos6)S 1 (iu, q)}. (26) 



V dVLduj 8tt 3 ' 

where E and E' are the incoming and outgoing (anti) neutrino energies, respectively. The 
transferred energy is u = E — E' and the transferred three-momentum is obtained from the 
(anti) neutrino incoming (k) and outgoing (/?) three-momentum as q = k — k' . The neutral 
current vector and axial vector charges are Cy = 1/2 and Ca = ~9a/^ where g a = 1.260 31|. 
Gf/(^c) 3 = 1.166 39(1) x 10~ 5 GeV~ 2 is the Fermi coupling constant. Once the response has 
been evaluated it is straightforward to evaluate the neutrino mean free paths in the medium 
using 

Let us mention that, although the low temperature case in an astrophysical scenario arises 
after the early cooling neutrino phase, it is worth sizing the dependence of the neutrino 
opacity on medium effects, spin polarization and strong magnetic fields with respect to 
the non-magnetized case in an analogous way to the study of the ideal cases of isospin 
symmetric nuclear matter and pure neutron matter evaluated before addressing the more 
realistic situation where different hadronic species in beta equilibrated matter arise. 



III. RESULTS 

In this section we present results arising from a consistent treatment [26] minimizing the 
Helmholtz free energy calculated in Hartree-Fock approximation under the presence of a 
strong magnetic field _Fm(A, p, B). We have considered magnetic fields up to a maximum 
internal value B max « 10 18 G according to the scalar virial theorem. In Fig. [TJ and Fig. [2] we 
plot the vector, S°(w,q), and spin, S" 1 ^, g), response functions respectively, as a function 
of the energy transfer, u, for the Skyrme SLy7 interaction at saturation density po m the 
RPA monopolar (solid line) and RPA dipolar (dashed line) approximations at a value of 
the transferred three momentum, q = 0.25 /m _1 . We can see that dipolar terms add non- 
negligible contributions to the monopolar response functions due to the fact that they retain 
the energy and momentum dependence of the residual interaction. In the spin density 
channel the effect is a shift in the collective mode to lower energies while in the vector 
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channel is the opposite. In Fig. [2] the spin response function shows the presence of spin 
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FIG. 1: Vector response function 5° at q = 0.25 fm~ l for a magnetic field B = 10 17 G at saturation 
density with the Skyrme SLy7 foce for I = and 1 = 1. 



zero sound modes as obtained in the dipolar approximation. Under these conditions there 
is a very mild polarization (A « —1%) of the neutron plasma. The collective spin zero 
sound modes are very close and their shift is hardly visible on the plot, uj + « oj_ « 33.0 
MeV. The ph limit for the cr-polarized population is given by u a = £(£±^£^1 j n pjg [2J we 
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FIG. 2: Spin response function S 1 for the same conditions as in Fig. [TJ 



can see the spin response function, S 1 , in the RPA dipolar approximation as a function of 
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the energy transfer, oj, for the Skyrme SLy7 interaction at saturation density, po, and zero 
temperature for different values of the transferred three momentum, q = 0.05 fmT 1 (solid 
line), q = 0.1 fmT 1 (short dashed line), q = 0.25 fm~ l (long dashed line) and q = 0.5 jm~ x 
(dotted line) respectively. For this plot we have set the magnetic field strength B = 10 17 
G. We can see that the ph continuum is achieved at a value max[uj + ,ujj\ and again at this 
density the very mild polarization of the system results in the relative separation of the 
modes to be less than 1% although increasingly larger as q increases. 
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FIG. 3: Spin response function S 1 in the RPA dipolar approximation for values (from left to right) 
of q = 0.05,0.1,0.25,0.5 fm~ l for a magnetic field B = 10 17 G at saturation density with the 



Skyrme SLy7 parametrization. 



In Fig. H] we can see the spin response function S* 1 as a function of the energy transfer, 
u>, for the Skyrme SLy7 interaction at p = 0.5p (upper pannel) and p = 3p (lower pannel) 
in the RPA monopolar (/ = 0) and RPA dipolar (1 = 1) approximations for a value of 
the transferred three momentum, q = 0.5 fmT 1 . In order to size the effect of a magnetic 
field we have studied limiting cases ranging from B = G to the maximum field strength 
considered in this work, B = 10 18 G. We can see that in the low density case there is a high 
energy mode shifted in the dipolar case with respect to the monopolar case. However the 
inclusion of magnetic field produces an undistiguishable (on the plot) shift in the strength 
due to the low induced polarization at this density 2(|. However, at higher density a low 
energy mode arises, since a ferromagnetic transition is near. The presence of magnetic field 
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FIG. 4: Spin response function at densities p = 
0.5po (upper panel) and p = 3po (lower panel) 
with the SLy7 model in the RPA Z = 0, 1 cases 
for limiting values of B. See text for details. 



tends to decrease (increase) the low (high) energy strength of the response slightly. At low 



densities a realistic treatment would include the possible existence o 
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25). 



magnetic pasta phases 



analogous to the clustered configurations obtained in [22, 

In order to further clarify the contribution of the magnetic field to the change in the 
spin response function we plot on Fig. [5] the ratio of change of the ph continuum for the 
cr-polarized population component with respecto to the B = 0G value, uj at a given value 
of q, that we define as R° = versus the logarithm (base 10) of the magnetic field 

strength for densities p/po = 1, 2, 3.4 in solid, long dashed and short dashed line respectively 
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and q = 0.25 frrT 1 in the RPA dipolar approximation. The change for fields below B pa 10 16 
G is almost negligible. We can see on the plot that the ratio is positive (negative) for R° for 
the parallel (antiparallel) spin populations in the p/p Q = 1,2 cases. In the p/p = 3.4 and 
due to the proximity of the ferromagnetic phase transition both ratios are negative and the 
antiparallel collective mode is largely affected to about ~ —11%. In Fig.EJand Fig.[7]we show 
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FIG. 5: Ratio of change of the ph continuum energy (in percentage) versus the logarithm of the 
magnetic field strength in the RPA dipolar approximation at different densities, p/po = 1,2,3.4 
and q = 0.25 fm~ l with the Skyrme SLy7 interaction. See text for details. 




the real and imaginary parts, respectively, of the dynamical Lindhard function, x^ 7 > as 
calculated with the Skryme SLy4 interaction at p = 3po with the RPA dipolar approximation 
for a value of B = 5 x 10 17 G and q = 0.1 fmT 1 . We can see that in the real and imaginary 
parts there is a shift to the lower energies in the dominant spin down component (~ 10%) 
and low energy modes are excited since at this density a ferromagnetic transition is near. 
For the imaginary part both the (+, +) and (— , — ) component are negative while the mixing 
part (+, — ) gives positive contribution so that the structure function in the S = 1 (S=0) 
channel is increased (decreased). As density grows the uj + and uj- modes drift away from 
each other since the spin polarization of the plasma is larger. 

In order to see the differences induced by the use of several parametrizations of the 
Skyrme force we plot in Fig. [8] the spin response function versus energy transfer for the 
Skyrme SLy4 (solid line) and SLy7 (dotted line) models at B = 5 x 10 17 G in the RPA 
dipolar approximation, p/po = 3 and q = 0.1 ]vrC x . There is a shift in strength and energies 
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FIG. 6: Real part of the Lindhard function versus energy transfer as calculated with the Skyrme 
SLy4 interaction with B = 5 x 10 17 G in the RPA dipolar approximation, p/po = 3 and q = 
0.1 fm,- 1 . 
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FIG. 7: Same as Fig. [6] but for the imaginary part. 



of different sign for both parametrizations in the high and low energy modes excited due 
to the fact that effective masses and polarization in the plasma are slightly modified for 



both interaction 



26j. Let us notice that the energy shifts belong to the sub-MeV range, 



the same order of the induced resolution of the degeneracy of the uj a modes in the polarized 
system. In Fig. [9] we have plot the ratio of the neutrino mean free paths Ri = Ab a ~ b A ^ =0 versus 
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FIG. 8: Spin response function versus energy transfer for the Skyrme SLy4 and SLy7 parametriza- 
tions at B = 5 x 10 17 G in the RPA dipolar approximation, p/po = 3 and q = 0.1 /to -1 . 

the logarithm (base 10) of the magnetic field strength in the RPA monopolar (I = 0) and 
dipolar (I < 1) approximation for the Skyrme SLy4 and SLy7 parametrizations at a density 
p/Po — 3 and and incoming neutrino energy of E u = 15 MeV. Since this ratio stays positive 
in the x-axis the neutrino opacity is smaller as magnetic field increases. We can see that 
for the Skyrme SLy4 parametrization at this density there is a negligible difference in the 
monopolar and dipolar case and the decrease in the opacity with respect to the B = case 
is ~ 4%. For the Skryme SLy7 parametrization we can see that the change in the dipolar 
case is smaller than in the monopolar case as magnetic field increases. At the maximum 
magnetic field strength considered at this density there is a 10% decrease in opacity with 
respecto to the non-magnetized case in the dipolar response. 

In order to check the density dependence of the relative changes induced by the inclusion 
of magnetic fields we plot in Fig. [10] the ratio of neutrino mean free paths in the RPA dipolar 
approximation computed for a fixed value of magnetic field stregth B = 10 17 G with respect 
to the B = G case, Rb = Xb=W \J 3 Xb= ° as a function of the density with the Skyrme SLy4 
(solid line) and SLy7 (dashed line) parametrizations for a value of the neutrino incoming 
energy E v = 15 MeV. We can see that the SLy4 predicts a larger decrease in the neutrino 
opacity than for the SLy7 case as the density grows to approach a ferromagnetic transition. 
This behaviour is in agreement with the dramatic decrease of the neutrino mean free path 
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FIG. 9: Ratio of neutrino mean free paths calculated with the SLy4 and SLy7 in the RPA approx- 
imation with I = 0, 1 as a function of the logarithm of magnetic field strength and E u = 15 MeV 
for p/po = 3. See text for details. 



as a phase transition arises in the system as obtained for the non-magnetized systems 
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FIG. 10: Ratio of neutrino mean free paths as a function of density calculated with the SLy4 and 
SLy7 parametrizations for a value of magnetic field strength B = 10 17 G and a neutrino energy 
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15 MeV. 



IV. SUMMARY AND CONCLUSIONS 



We have investigated in the context of the Landau Theory of normal Fermi Liquids, 
the effect of a strong magnetic field on the density (spin density) fluctuations in a two- 
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component partially magnetized pure neutron system within the framework of the non- 
relativistic Hartree-Fock approximation. To describe the nuclear interacting system we 
have used Skryme effective forces under the form of the SLy4 and SLy7 parametrizations. 
Due to the tiny value of the neutron magnetic moment we have considered magnetic field 
strengths up to the limiting value B pa 10 18 G as allowed by the scalar virial theorem. The 
matter vector and vector-axial response functions, S , ^ s=0 ' 1 ^ are obtained from the imaginary 
part of the Lindhard function, calculated in the RPA approximation using monopolar (I = 
0) and dipolar (I < 1) expansion terms in the polarized quasiparticle interaction matrix 
element. The dipolar approximation is exact for the Skryme force and adds non-negligible 
contributions to the usually considered monopolar RPA structure functions, shifting the 
plasma collective modes to the lower (higher) energy range of the spectrum in the vector 
(vector-axial) functions with respect to the monopolar case. The collective modes, u a , 
depend on the net polarization so that when it is not vanishing the degeneracy of the 
modes in the non-magnetized case is resolved into pairs uj + ,uj^. When the polarization is 
mild there is a sub-MeV energy mode separation increasing as density and magnetic field 
strength grows. As most of Skyrme models, near the predicted ferromagnetic transition a 
low energy mode arises, signaling the energetic instability towards the A = — 1 case. 

In order to size the contribution to the results of the inclusion of a magnetic field we have 
also computed the ratios of variation of the collective modes with respect to the B = G 
case finding that the as density and magnetic field strength increases, the relative splitting 
grows. Using different Skyrme parametrizations in the NN interaction produces a sub-MeV 
shift in the collective mode location. Finally, to understand the response to a weak probe 
interacting via neutral currents we find that the magnetic field tends to decrease neutrino 
opacities giving a non-negligible ~ 10% correction or higher for a density of about three 
times saturation density at the maximum magnetic field strength log 10 B(G) = 18 and 
incoming neutrino energy E v = 15 MeV. The inclusion of finite temperature magnetized 
matter in the presence of magnetic fields should be taken into account in astrophysical 
scenarios by considering beta-equilibrated matter so that polarization of protons would 
include corrections to the present result. This task should be addressed in future work. 
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